Specific heat at constant volume in the thermodynamic model 
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A thermodynamic model for multifragmentation which is frequently used appears to give very 
different values for specific heat at constant volume depending upon whether canonical or grand 
canonical ensemble is used. The cause for this discrepancy is analysed. 
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The motivation for this work is just one puzzle. A thermodynamic model, often used for fitting data, appears to 
give very different answers for specific heat at constant volume depending upon whether the canonical or the grand 
canonical ensemble is used. We wish to resolve this issue. The relevant papers are [1] and [2]. Although for practical 
applications a much more sophisticated two component version of the model is used, [3-5] here, as in [1,2] we use 
[ one kind of particle as in the original formulation [1]. If a system has A nucleons in a volume V at temperature 
^sj ■ T the system can break up into various composites which have a binding energy consisting of volume and surface 
energies. Excited states of the composite can also be included. Assuming that the break up takes place only according 
to availability of phase space, it was shown in [1] that this problem for A nucleons can be solved numerically with 
arbitrary accuracy very easily. One can also solve the problem in a grand canonical ensemble [2,6]. We do not give 
any details here as they are given in many places including [1,2]. 

By specific heat we will always mean specific heat per particle. It was shown in the original paper [1] that for a 
fixed A, the total number of particles and a fixed V , the specific heat as a function of temperature went through 
a maximum at a certain temperature which was labelled the boiling temperature. A numerical calculation which is 
rather easy and can be made sufficiently accurate showed that for a fixed p = A/V if we increase the number of 
' particles the height of the peak rises and the width decreases. Since one sees no reason why this behaviour should 
' change at certain high value of A, it was concluded that the specific heat behaves like a delta function in the limit 
. A — > oo. There is a physics picture one can associate with this. Below the boiling temparature there is a blob of liquid. 
' Just above the boiling temperature this blob changes into a system of smaller composites and nucleons. Qualitatively 
a delta function would emerge if a finite fraction rather than an infinitesimal fraction of the blob converts into gas 
with an infinitesimal increase of temperature. 

In Fig. 1 we show the specific heats in the canonical model for p — /3o/2.7 for A — 200 and A ~ 2000. All calculations 
shown here will be at this p. The grand canonical results are also shown in the figure. For grand canonical we solve 
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where km is the number of nucleons in the largest cluster allowed in the system. Here 
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To do the grand canonical calculation for A=2000 we set km in the above equation at 2000. The average value of 
< nfe > is then given by, 

< rtfc >= exp(fc/?/i)(Dfc X V (3) 

where V is the appropriate value of 2000 nucleons at freeze-out, i.e., V = 2000 x 2.7/pq. Denoting < > as the 
average number of composites which has k nucleons (i.e., the average number of trimers is < ns >) we have 

2000 

J2k<nk>= 2000. (4) 

k=l 

It should be realised that all A's (upto oo) are included in the grand canonical ensemble but setting fcm=2000 in eq. 
(1) signifies that the largest cluster has 2000 nucleons. In terms of canonical partition functions one has 
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Qgr.can = ^ ex.p{(3lJ.K)Q K,km (5) 
K=0 



where QK.km i'' the canonical partition function of K nuclcons but with the restriction that the largest cluster has 
only km{= 2000) nucleons. The quantity /3/x is known from solving eq.(l) with ^^=2000. Similar arguments hold 
when A;to=200 but we will not discuss this case and concentrate only on A = km = 2000. 

Fig.l shows that the specific heats in the canonical and grand canonical models are very different. In both the 
models, the peak value of the specific heat increases when we go from 200 to 2000 particles and the widths decrease 
but the results are much more dramatic in the canonical model. Since A=2000 is a large number in the context of 
nuclear physics, we try to understand the cause of this difference. In particular it is not obvious that the specific heat 
in the grand canonical model will attain extraordinary heights and/or miniscule widths. We will show that the cause 
of discrepancy between the canonical and grand canonical models is the very large fluctuation in the particle number 
in the grand canonical ensemble. We can investigate this in two ways, one more detailed than the other. We note 
that 
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and we can calculate fluctuation exploiting the well-known relation 
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= ^ < rife > (7) 
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But we can also exploit the fact that we know Q/f,fc„ upto rather large values of K and the value of /J/x from the 
grand canonical calculation. Hence we can use the equation (5) also to calculate fluctuation. For practical reasons, 
the upper limit of K will have to be cut off. The upper limit of K in the sum above was 10,000. Since we are 
investigating A=2000 one might a priori assume this should be adequate. 

The fluctuations calculated from eqs.(5) and (7) are shown in Fig. 2. One sees that there is a temperature above 
which the fluctuations are small. At these temperatures the grand canonical value of speciflc heat is indistinguishable 
from the canonical value. But as temperature is lowered, fluctuations grow rapidly and the results begin to diverge. 

It is interesting to study fluctuations further. The probability of K particles being in the grand canonical ensemble 
is oc e^/^M+inQK ^Eq_ 5) g^^^j .^g pjg^ pjg 3 exp[/3/i(if - A) + IhQk - ItiQa]- This takes the value 1 at ii" = A and 
in the normal picture of grand canonical ensemble, would drop off rapidly on either side of A. This does happen at 
temperature higher than the boiling temperature. The case at temperature T=7.7 McV corresponds to a standard 
scenario. But the situation at temperature 7.3 MeV is drastically different. The probability does not maximise at 
K = A but at a lower value. It is also very spread out with a periodic structure. At temperature 2.0 MeV, the 
probability of having no particle is higher than the probability of having K = A. We notice that here also there is a 
periodicity in the probability distribution. The periodicity is 2000 and is linked with the fact that in the case studied 
the largest composite has 2000 nucleons and at low temperatures, this composite will play a significant role. 

We can now understand why the speciflc heat curve is so flat in the grand canonical ensemble in Fig 1. Even though 
the average number A is 2000, the ensemble contains large components of K < A (thus have lower density and peaks 
of specific heat below 7.1 MeV) and K > A (which have peaks at higher than 7.1 MeV). It is this smearing which 
makes the specific heat peak much lower and much wider. 

In Fig. 4 we have shown canonical and grand canonical results for the total energy of 2000 particles. The canonical 
result suggests that starting from a low temperature, energy increases at a finite rate (implies a finite value of Cy), 
followed by a sudden rise (will lead to infinite specific heat), followed again by a regular behaviour. Thus the transition 
region is marked by two different values of Cy with a delta function switched in between. In the grand canonical 
model for 2000 particles one sees only the discontinuity in the value of Cy. If one was calculating the Cy directly in 
either the canonical or grand canonical model and for a very large system there is indeed a delta function in Cy , one 
will miss the delta function since it has zero width. To see that there is indeed one, we should instead calculate the 
total energy and see that at a given temperature T there is a huge difference in the total energy for an infinitesimal 
increase in T. 

One might get the impression that for finite systems the use of a grand canonical ensemble is very dangerous. For 
many observables it is quite acceptable. However, a recent work shows [7] that interpreting data according to grand 
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canonical ensemble can lead to a serious error in estimating temperature. It thus depends on the particular observable 
being calculated. 

Returning briefly to the interesting periodicity seen in Fig. 3, it arises because at low temperatures it is advantageous 
for the system to form large clusters. In the example studied in Fig. 3 where the largest cluster size was 2000, at 
temperature 2 MeV when the system has 8000 nucleons we have essentially 4 clusters, each of size 2000 (we are 
confining ourselves to canonical calculation of course). This becomes less precise at higher temperature. For example 
at 6 MeV temperature if the system has 4100 nucleons nearly 4000 of them are distributed in large clusters (more 
precisely, on the average, 3950.22 nucleons are bound in clusters of sizes between 1800 and 2000) and nearly 100 
of them are in lighter clusters (149.78 on the average). If we now go to a system of 6100 nucleons, on the average 
5940.22 are bound in very heavy clusters and 159.78 in lighter clusters. Approximately speaking, at low temperatures, 

ZnQx+2000 ~ InQK + ln{Vui2ooo)- 

The low temperature periodic structure at T = 2 MeV can be qualitatively understood using the following resiilts. 
As T ^ 0, the system will go to the largest cluster allowed, and in this case, km = 2000. For example, a.t K = 10, 000 
(the total number of nucleons in the system) a result of 5 clusters of size 2000 follows. The mean number of cluster 
of size km{= 2000) is: 



Also the factorial moment < (nfe^ — 1) > is given by. 
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And in general. 



< nfe„(nfe„ - l)...(nfe„ - n + 1) > = (wfc„)" ( ^^^^^ ) . (8) 



Thus the points at 8000, 6000, 4000, 2000 and can be related to the factorial moments of the Uk^ distributions as 
T ^ 0. At T = these are 5x4x3x2x1, 5x4x3x2, 5x4x3, 5x4 and 5 at iC = 0, 2000, 4000, 6000 and at 
8000 respectively. The heights of the peaks in Fig. 3 will be determined by, 

^ < nk^iuk^ - l)...(nfe„ - n + 1) > ^e^^'(^-^^+^)^^) 

as T — > 0. Even at T = 2 MeV, this expression is very accurate. 
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FIG. 2. Fluctuations calculated using eqs.(5) and (7). The solid line corresponds to using eq. (5) with K cut off at 10,000 
and the dotted line corresponds to using eq. (7). Ti correspond to the temperature where the specific heat maximises in the 
canonical calculation and T2 to the temperature of highest specific heat in the grand canonical calculation. 
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FIG. 4. Caloric curves in the canonical and grandcanonical models, for a system of 2000 particles. 
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